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Abstract 

Many challenging Graver bases computations, like for multi-way tables in statistics, have a 
highly symmetric problem structure that is not exploited so far computationally. In this paper we 
present a Graver basis algorithm for sublattices of Z n that exploits existing symmetry. 

1 Introduction 

Graver bases, originally introduced by Graver |3] for use in integer programming, have a variety of 
interesting applications. Besides providing improving directions for integer programs f° r 
stochastic integer programs UJ, and even for certain convex integer programs 0E31; Graver bases can 
for example also be used as Markov bases for sampling in statistics or as a superset from which a 
universal Grobner basis of the toric ideal I a '■= (% u — x v : Au = Av, u,v € Z" } can be extracted \1'2\ . 

Unfortunately, the size of Graver bases increases quickly with the dimension, making it very hard if not 
impossible to compute them in practice. In several applications, however, as for example in algebraic 
statistics, the problems involve a high symmetry that should make it much easier to compute the 
Graver basis in terms of (relatively few) representatives. In HHH], the authors exploit existing symmetry 
of lattices that arise from a single constraint of the form a T x = (mod p), p € Z + , in a similar way 
as we do in Lemma Tl. 21 and in Corollary 12 .41 below. 

Let us start our presentation by defining the notion of a Graver basis and by giving an example that 
demonstrates the problem we are interested in. 

The Graver basis <7(A) associated to a lattice A C Z™ consists exactly of all □-minimal nonzero 
elements in A, where for u,v £ Z™ we say that u C v if u^v W > and < \v^\ for all 

components j = 1, . . . ,n, that is, if u belongs to the same orthant as v and its components are not 
greater in absolute value than the corresponding components of v. Note that Graver originally defined 
this set only for the case A = ker(A) C\Z n for given matrix A 6 Z dx ™, but his definition can be readily 
extended to the definition we gave above. 
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Example 1.1 Consider the set of all 3 x 3 tables/arrays whose entries are filled with integer numbers 
in such a way that the sums along each row and along each column are 0. One particular example is 
the table 



If we encode the 9 entries of the table as zi, . . . , zg, then the set of 3 x 3 tables coincides with the 
integer vectors in the kernel of the matrix 
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that is, with all z G Z 9 satisfying A^x^z = 0. As defined above, the Graver basis of ^3x3 consists of 
all C-minimal nonzero tables among them. The particular 3x3 table above does not belong to the 
Graver basis of ^3x3, since 



Using the computer program 4ti2 0], we find that the following 15 vectors (and their negatives) 
constitute the Graver basis of ^3x3: 
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However, there is an obvious symmetry group S3 x S3 x S2 operating on the set of 3 x 3 tables 
whose elements transform a given table v G ker(A3 X 3) into another table w G ker(A3 X 3) by suitably 
rearranging components (permuting rows or columns, flipping the table along the main diagonals). If 
we take these symmetries into account, we see that among these 15 elements there are in fact only 
two essentially different elements: 

(1, -1, 0, -1, 1, 0, 0, 0, 0) 
(1, -1, 0, -1, 0, 1, 0, 1, -1) 

or, in a more array-like notation: 

/ 1 -1 \ / 1 -1 \ 

-1 1 and -1 1 . 
v / \ 11/ 

It should be clear that for bigger or for higher-dimensional tables, this difference in sizes becomes far 
more striking, since the acting symmetry groups are much bigger. □ 



In the following, let a sublattice A C Z™ be given and let 5a E S n be a group of symmetries such that 
for all v G A and for all a G Sa we have that also a(v) := (v"^, . . . ,v a ^) G A. Finally, denote by 
orbs A (w) := {cr(v) : a G Sa} the orbit of v under 5a. 

Lemma 1.2 Let A C Z™ be a lattice, let Sa E S n be the group of its symmetries, and letg,g',s,s' G A. 
Then the following holds: 

• If 9 E s then a(g) C a(s) for every a G 5a. 

• If there is no g' G orbs A (g) with g' Q s, then for every s' G orbs A (s) there is no g" G oibs A (g) 
with g" C s'. 

• If v (A) then orb SA (w) C g(A). 

Proof. The first statement follows immediately from the definition of C. 

For the second statement, let s' = a(s) for some a G 5a and assume there is some g" G orbs A (<?) 
with g" C s'. Then g' := o-^ 1 (g") C cr _1 (.s') = s and g' G orbs A (g). A contradiction to the assumed 
non-existence of such g' . 

For the last statement we have to show that with v G 5(A) also the full orbit orbs A (w) lies in 5(A). 
For this it suffices to assume that there is some a(v) G orbs A (v) that could be written non-trivially 
as a(v) — Wi + u>2 with W\,VJ2 E cr(v). But this would imply v — a^ 1 (w\) + cr~ 1 (w2) with nonzero 
o-~ 1 {w 2 ) E v, which contradicts v G 5(A). Thus, a(v) must belong to 5(A). □ 

As a consequence of this lemma, 5(A) decomposes completely into full orbits. Our task of computing 
5(A) thus reduces to computing representatives of these orbits, and to collect them into a set 5 sym (A). 
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By the previous lemma, we recover £7 (A) via 

0(A) = |J orb s ». 

Note that this last expression does not compute a superset of 5(A). In contrast to this, the last 
statement of Lemma 11.21 fails to be true in general for minimal toric Grobner bases or for minimal 
Markov bases associated with the lattice A. 

2 Computing the Graver Basis 

In this section, we adapt Pottier's algorithm |11| to deal with the symmetries of A in the computation 
of 5(A). Note that this algorithm is not the fastest way to compute Graver bases directly. However, 
by adapting this algorithm, it will be easier for us to exploit the given symmetries (and to present the 
main ideas). The state-of-the-art algorithm that is based on the positive sum property of Graver bases 
|Sj needs to break the symmetry, see Section [3] Nonetheless, we show how to exploit the symmetries 
also in this situation and arrive at an even faster "symmetric" algorithm. 

Algorithm 2.1 (Algorithm to Compute G(A)) 
Input: generating set F of A over Z 
Output: a set G which contains G(A) 

G := F (J —F 

C:= U {/ + <?} 

f.gec 

while C 7^ do 

s := an element in C 

C:=C\{s} 

f := normalForm(s, G) 

iff ^0 then 

G:=GU{f} 

c-.= cu U {f + g} 

gee 

return G. 



Behind the function normalForm(s, G) there is the following algorithm. It aims at finding a relation 
s = gi + . . . + g r , gi £ G, Qi C s, i = 1, . . . , r. The function normalForm(s, G) returns if it found 
such vectors £ G, or it returns a vector / £ A such that such a desired relation exists with elements 
from GU {/}. 
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Algorithm 2.2 (Normal Form Algorithm) 

Input: a vector s, a set G of vectors 

Output: a normal form of s with respect to G 

while there is some g £ G such that g C s do 

s := s-g 
return s 

This algorithm due to Pottier always terminates and the set 5(A) is exactly the set of all C-minimal 
vectors in the final output G. Termination is guaranteed by the following lemma, which we will employ 
again later. 

Lemma 2.3 (Gordan-Dickson Lemma, Sequence version, JTjj) 

Let {p\,p2, ■ ■ ■} be a sequence of points in Z™ such that pi pj whenever i < j. Then this sequence is 
finite. 

Now let us adapt the algorithm to exploit the symmetries. Let us start with an immediate consequence 
of Lemma fl .21 

Corollary 2.4 If s = g x + . . . + g r , gi C s, i = 1, . . . , r, then a(s) = a(gi) + . . . + a(g r ), a(g t ) C a(s), 
i = 1, ...,r for every a S 5a- Thus, if G = [J ovbs A (v) and if normalForm(s, G) — then 

normalForm(cr(s), G) = for every a € 5a- 

In other words, if a representation s = g± + . . . + g r , 9i 6 G, gi C s, i = 1, . . . , r has been found, the 
symmetry of A already guarantees existence of a similar representation for every element in orbg A (s). 

Finally, we are in the position to exploit symmetries in Pottier's algorithm. The main differences to 
the original algorithm will be that instead of keeping the sets G and C in memory, we only store 
their representatives under the given symmetry in sets G sym and C sym . (At any point during the 
"symmetric" algorithm we may go back to the original algorithm by replacing all elements in G sym 
and C sym by the vectors from their orbits under iSa-) 

Moreover, there a few more changes. Once we have found a nonzero vector / that is to be added to 
G (and to G sym as a new representative), we assume that we add the full orbit orbs A (/) to G, as 
we know that the Graver basis would contain the full orbit if / was in fact a Graver basis element. 
Accordingly, instead of adding only the vectors 

gee 
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to G, we immediately include the vectors 

U {/'+£}■ 

/'Gorbs A (/),geG 

As G will always be a union of full orbits, this last expression can be transformed to 

(J {f + 9}= (J U if + 9'}= U U OTb s A (.f + £')• 

/'eorb SA (/), 9 eG 9GGaym /' e orb SA (/) seG^" 9 'eorb SA ( 9 ) 

g' e orb SA (3) 

Therefore, we update G sym as follows: 

c ,sym = c ,sym u (J (J + 

9 £G"»j'arb Sl ( S ) 

Note that since 

|J orb SA (/ + g')= |J orb SA (.f+. 9 ) 

g'eorbs A (g) /'£orbs A (/) 

we have a choice in adding either all vectors {/ + <?'} or all vectors {f' + g} to G sym . Clearly, one would 
choose to add as few new representatives to G sym as possible to keep the number of S-vectors that 
need to be reduced small. After all, these reductions are the most expensive part of the algorithm. In 
the following algorithm, rcps A (H) for a set H C A of vectors shall denote a set of representatives of 
H under the symmetry group 5a. 

Algorithm 2.5 (Algorithm to Compute G sym (K) ) 
Input: generating set F of A over Z 
Output: a set G which contains 5(A) 

Qsym ._ rcpg A (_F U —F) 

C^:=rep SA ( U {/ + <?}) 
\i\geG ) 

while C sym 7^ do 

s := an element in C sym 

Qsym ._ (jsym \ r g \ 

f := normalForm(s, G sym ) := normalForm(s, G) 
iff^O then 

Qsyrn ._ Qsym U {/} 

(jsym ._ (jsym y (J (J |j + g /| 

9eG 3 'eorbs A (9) 

refarn G = orbs A (G 52/m ) . 



G:=orb SA (F U-F) 
C:= U {/ + <?} 



G:=G\orb SA (s) 
/ := normalForm(s, G) 

G:=GUorb SA (.f) 
G:=GU U {f' + g} 

/'Sorbs A (/), 9 £G 
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In this algorithm, wc compute normalForm(s, G sym ) via normalForm(s, G sym ) := normalForm(s, G). 
Clearly, from a practical perspective, one would not want to keep the huge set G in memory. Then, of 
course, one needs to think about how to compute normalForm(s, G sym ) efficiently if only G sym instead 
of G is available. (For example, G might be simply too big to be kept in memory.) This is still an open 
question and any significant improvement in the solution of this problem would lead to an equally 
significant improvement of the overall algorithm. 

Lemma 2.6 Alaorithm \2.5\ always terminates and returns a set G containing Q(A). 

Proof. To prove termination, consider the sequence of vectors in G sym \reps A (-FU — F ) = {/i, /2, . . .} 
as they are added to G sym during the run of the algorithm. By construction, we have fi % fj, that 
is (/+, /r) ^ (fj, fT) whenever i < j. Thus, by the Gordan-Dickson Lemma, this sequence must be 
finite and the algorithm terminates. 

Note that throughout the run of the algorithm, we always have G = orbs A (G sym ). Upon termination 
we know that normalForm(/ + g, G) = for every pair of vectors f,gE& Thus, G must contain the 
Graver basis Q(A). □ 



3 Computing the Graver basis faster 

In this section we introduce a special generating set of A. This set not only decreases the number 
of sums f + g that need to be added to G, but more importantly, it reduces the amount of work 
to compute normalForm(s, G) tremendously. As the latter computation is the most expensive part 
of Pottier's algorithm, this heavily speeds up the computation of <?(A). Moreover, we introduce a 
so-called critical-pair selection strategy that chooses the next element s from G according to a certain 
rule. This, together with the special input set, will imply that the set G returned by our algorithm is 
exactly the Graver basis £(A). Not a single unnecessary vector is computed! 

In Pottier's algorithm, one has to wait until the very end to extract the Graver basis from the returned 
set G. In contrast to this, our algorithm provides a certificate of C-minimality for each vector that is 
added to G. Consequently, at any point during the computation, a subset G of G(A) is known. 

Let us assume from now on that A is generated by d vectors and that the first d components of 
the vectors of A are linearly independent, that is, the only d-dimcnsional vector orthogonal to the 
projection of A to the first d components is the zero vector. (This condition can easily be achieved 
algorithmically on a lattice basis of A by integer row operations and by switching components.) Let 
7r denote the projection of an n-dimensional vector onto its first d components. As the last n — d 
components of A are linearly dependent on the first d components, a vector tt(v) 6 7r(A) can be 
uniquely lifted back to v € A. Moreover, this can easily be done algorithmically. 

Next let us define a norm |.| on vectors v S A as follows: ||u|| :— ||7r(u)||i, where ||.||i denotes the 
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Li-norm on R d . It can easily be checked that, under our assumptions on A, this defines indeed a norm 
on A. It is this norm definition that breaks existing symmetry of the given problem. 

Finally, we are ready to state the algorithm behind the implementation in 4ti2 that seemingly defines 
the current state-of-the-art in the computation of Graver bases. 

Let 7r(F) denote set of all □-minimal nonzero vectors in 7r(A). Note that n(F) has the positive sum 
property with respect to 7r(A), that is, every vector n(v) E 7r(A) can be written as a positive integer 
linear combination tt(v) = a i n {fi) of vectors 7r(/j) from the set n(F) that lie all in the same orthant 
as ir(v). In other words, on 6 Z >0 and 7r(/j) C tt(v) for all i in this linear combination. Clearly, this 
nice integer linear combination n(v) — ^a^ifi) can be uniquely lifted to v = ^ctifi showing that 
F is in particular a generating set of A over Z. 

Moreover, this set n(F) can be computed for example via Pottier's algorithm once a lattice basis F 
of 7r(A) C Z d over Z is given as input. Note that all these C-minimal vectors in n(F) must lift to 
C-minimal elements in A, since they are already indecomposable/minimal on the first d components. 
Therefore, the computation of tt(F) is usually far less expensive than computing <?(A) itself. In our 
computational experiments, it often happened that this set n(F) simply contained the unit vectors in 
Z d and their negatives. 

Algorithm 3.1 (Faster Algorithm to Compute G(&)) 

Input: set F C A such that n(F) is the set of ^-minimal nonzero vectors in 7r(A) 
Output: a set G which contains Q(A) 

G:=F 

C:= U U + 9) 

/,<?GG:7r(/) and Tr(g) lie in the same orthant of R d 

while C 7^ do 

s := an element in C with smallest ||.|| -norm 

C:=C\{s} 

f := normalForm(s, G) 

iff ^ then 

G:=GU{/} 

C:=CU U {/ + .?} 

pGG:7r(/) and ir(g) lie in the same orthant of R d 

return G. 



Due to our special input set, the function normalForm(s, G) can be sped up as follows. 
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Algorithm 3.2 (Faster Normal Form Algorithm) 

Input: a vector s, a set G of vectors 

Output: a normal form of s with respect to G 

if there is some g £ G such that g C. s return 
return s 

Once a vector g £ G with g C s is found, a representation s = gi + . . . + g r , gi £ G, <?; C s, 
i = 1, . . . , r, must exist. Therefore, the function normalForm(s, G) can return immediately without 
explicitly constructing such a relation. 

Since we started with a very special input set, only those pairs of vectors f,g £ G lead to a critical 
vector in G, for which the projections ir(f) and 7r(g) lie in the same orthant of M. d . Finally, let us 
prove our claims. 

Lemma 3.3 Alaorithm \H.l\ always terminates and returns a set G containing G(A). 

Proof. To prove termination, consider the sequence of vectors in G\F = {fx, fi-, ■ ■ ■ } as they are added 
to G during the run of the algorithm. By construction, we have fi % fj, that is (f^, /r) ^ (/• , /~) 
whenever i < j. Thus, by the Gordan-Dickson Lemma, this sequence must be finite and the algorithm 
terminates. 

To prove correctness, let us assume that z £ Q(A) is not contained in the final output set G of 
Algorithm V3. II Without loss of generality we may assume that z has a smallest norm ||z|| among all 
such vectors from 5(A). Therefore, we can assume that all Graver basis elements g with ||<7|| < ||z|| are 
contained in G. In the following, we construct a contradiction to the assumption z g" G and correctness 
of Algorithm 13. H is proved. 

As F is contained in G, there is a representation z = ^ aiVi with positive integers a% and vectors 
Vi £ G with Tr(vi) C k{z). From the set of all such linear integer combinations choose one such that 
^2 oii 1 1 Vi 1 1 1 is minimal . 

Let us assume first that X) a ill u i||i > ll z lli- Therefore, there have to exist vectors u^,!^ in this 
representation which have some component k = fco of different signs. By construction, fco > d, as the 
Vi have all the same sign as z on the first d components. 

The vector Vi 1 + Vi 2 was added to C during the run of Algorithm 13.11 (as ^(u^) and tt(vi 2 ) lie in 
the same orthant of M. d by construction). If \\vi x + Vi 2 \\ = \\z\\, then all other Vi, i ^ ii,«2) must 
satisfy ||i>j|| = as ir(vi) □ 7r(z) for all i. But ir(vi) — implies Vi — and thus v^ + Vi 2 = z. 
Since v^ + Vi 2 (— z) is a vector that was added to G during the run of the algorithm, the vector z is 
eventually chosen as s £ G. Being □-minimal, we have normalForm(z, G) — z and thus, z must have 
been added to G, in contradiction to our assumption z ^ G. 
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Therefore, we may assume that \\vi 1 + Vi 2 \\ < \\z\\. However, since all Graver basis elements v with 
norm ||u|| < ||z|| are assumed to be in G, there must exist a representation + vi 2 = Pj v j f° r 
finitely many /3j £ Z>o, v'j £ G, and (3jVj C + v- l2 for all j. This implies that we have for each 
component k = 1, . . . , n, 

E« w i = IE Vi w i - ik < + i<i fe) i, 

3 3 

where the last inequality is strict for k = kg by construction. Summing up over k — 1, . . . , n, yields 
S j^jII^j" 111 = \\ v ii + v i2 Hi < IkJi + ll ll i 2 Hi- But now z can be represented as 

z = a^v^ + a l2 v l2 + E 

= ^ Pjv'j + (a h - l)v h + (a i2 - l)v i2 + E 

and it holds 

E&KHi + On -i)Kiii + -i)Kih+ E ^n^iii <E a *Hi i 

in contradiction to the minimality required on ^ ci!i]]i>j||i- Thus, our assumption X) a «ll w illi > ll z lli 
was wrong and ^ill^illi — IMIi must hold. 

But this last equation implies that Vi Q z for all i, contradicting C-minimality of z unless the rep- 
resentation z = \^2, a i v i is trivial, that is z = v± £ G. This, however, again contradicts our initial 
assumption z g 1 G and thus 5(A) C G. □ 

It should be noted that we did not make use of our selection strategy to prove termination and 
correctness of Algorithm 13.11 It is the following Lemma that provides a certificate for □-minimality 
of vectors in G. 

Lemma 3.4 The set G returned by Alaorithm \3.1\ equals 5(A). 

Proof. The main observation needed in this proof is that the norms ||s|| of the vectors s that are 
chosen from G form a non-decreasing sequence. This follows from the definition that normalForm(s, g) 
only returns either or s and from our condition that only vectors / + g are added to G whose first 
d components have the same sign pattern. The latter implies ||/ + .g|| = ||/|| + ||<?||. 

Now assume that some z £ G is not contained in 5(A). Thus, there is some g £ 5(A) with g C z. 
Under our assumptions on A, g C z implies \\g\\ < \\s\\. Since G contains 5(A) (and thus in particular 
the vector g) at the end of the algorithm and since after z only vectors / are added to G that have 
a norm ||/|| > ||z||, the vector g must have been contained in G already at the time when z was 
added to G. This however, implies that Algorithm \'6 . II must have computed normalForm(z, G) = 0, a 
contradiction to the assumption that z was added to G. □ 
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4 Computing the symmetric Graver basis faster 



Although the definition of |.| in the previous section broke most if not all existing symmetry in the 
problem, we will now combine the ideas of Sections[21and[!5]to a faster algorithm to compute symmetric 
Graver bases. The main idea is to use the norm |.| defined on Z" to define a norm on orbits: For 
T C A, we define ||T|| := min{||u|| : v € T}. Then the new symmetric Graver basis algorithm looks as 
follows. 

Algorithm 4.1 (Faster Algorithm to Compute G sym (A)) 

Input: set F C A such that ir(F) is the set of ^-minimal nonzero vectors in 7r(A) 
Output: a set G which contains G(A) 

Qsym ._ r ep SA (F) 

CW» : =re PSA ( U if + 9}) 
\f,geG ) 

while C sym 7^ do 

s := an element in C sym with smallest value of || orbs A (s)|| 

Qsym ._ Qsym \ r g \ 

f := normalForm(s, G sym ) 
iff^O then 

Qsym ._ Qsym \J ff\ 

C sym ._ C sym u y y {f + g'} 
g£G 9'Gorbs A (g) 

return G = orb SA {G sym ). 

Due to our special input set and our norm defined on orbits, we can again simplify and speed up 
the normal form computation normalForm(s, G sym ) := normalForm(s, orbs A (G sym )) by using Algo- 
rithm |^| instead of Algorithm 12.21 Again, from a practical perspective, one would want to compute 
normalForm(s, G sym ) without recovering or storing the huge set orbs A (G sym ). 

On the other hand, it should be noted that we do not, as in Algorithm 13. II have an orthant condition 
in Algorithm 14.11 that reduces the number of vectors added to C sym . This is done to "undo" the 
symmetry breaking caused by ||.||. 

Lemma 4.2 Algorithm \4-l\ always terminates and returns the set G = <?(A). 
Proof. The subsequent proof follows similar lines as the proof of Lemma 31 
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To prove termination, we consider again the sequence of vectors in G sym \ reps A (-F') = ft, ■ ■ •} as 
they are added to G sym during the run of the algorithm. By construction, we have fi % fj, that is 
(fi~'fr) ^ (fj>fj~) whenever i < j. Thus, by the Gordan-Dickson Lemma, this sequence must be 
finite and the algorithm terminates. 

Next we show that 5(A) C G. Assume on the contrary that this is not the case. Among all elements 
z € 5(A) \G choose one with || orbs A (z)|| smallest. Moreover, we may assume that z is a representative 
of orbs A (.z) with ||7r(z)||i = || orbs A (z)||. By our generating assumption on the set F C G, there exists a 
non-trivial representation z = Y] aiVi with positive integers a, and vectors Vi 6 G with ir(vi) C n(z). 
From the set of all such linear integer combinations choose one such that X^zll^lli 1S minimal. 
Note that from Tr(vi) C ir(z) and the fact that the relation z = Y2 a i v i i s non-trivial, we conclude 
IK(«;)||l < IK(z)||i and thus ||orbs A (u;)|| < ||orb SA (z)|| for alii 

Let us assume first that X] a ill l 'i||i > ll 2 lli- Therefore, there have to exist vectors Vi x ,Vi 2 in this 
representation which have some component k = fco of different signs. By construction, fco > d, as the 
Vi have all the same sign as z on the first d components. 

The orbit of the vector + Vi 2 was added to G sym during the run of the algorithm. If \\vi 1 + Vi 2 1| = ||z||, 
then all other Vi, i ^ i\, %2 must satisfy \\vi\\ = as 7r(«j) C ir(z) for all i. But n(vi) = implies Vi = 
and thus +Vi 2 = z. Since +Vi 2 (= z) is a vector whose orbit (or better: a representative of it) was 
added to C sym during the run of the algorithm, a representative of the orbit of the vector z is eventually 
chosen as s 6 G sym . Being □-minimal, we have normalForm(reps A (orbs A (^)), G sym ) = reps A (orbs A (z)) 
and thus, reps A (orbs A (.z)) must have been added to G sym , in contradiction to our assumption z ^ G. 

Therefore, we may assume that \\vi 1 +Vi 2 1| < ||z||. However, since all Graver basis elements with norm 
strictly smaller than ||z|| are assumed to be in G, we may continue literally as in the proof of Lemma 
rewrite z = Y] a^Vi and arrive at a contradiction to the minimality of ^ OjWj. 

Thus, we must have X) a ill w i||i = ll z lli' which implies that m C z for all i. As z is C-minimal, this is 
only possible if the relation z = Y] otiVi is trivial, that is z = V\. As v\ € G we now also have z G G 
as desired. 

To show that not only Q(h) C G but in fact 5(A) = G is true, assume G \ <7(A) ^ 0. Among all 
elements z £ G\Q(K) choose one with || orbs A (^)|| smallest. Moreover, we may again assume that z is 
a representative of orbs A (z) with ||7r(z)||i = || orbs A (^)||. Note that by Lemma IT7S1 and since z $ 5(A), 
not a single element in orbs A (.z) belongs to 5(A). 

We now start with a non-trivial representation z = Yl a^i with positive integers a.i and vectors Vi 6 
F C G with 7r(«j) C 7r(z). Clearly, 7r(u,;) C 7r(z) implies 1 1 vr (ui ) 1 1 1 < ll 7r ( z )lli an d thus || orbs A (^i) || < 
|| orbg A (z)|| for all i. Next, we follow similar steps as above (or more precisely: as in the proof of 
Lemma l3.3|) to change z = Y] onvi into a representation z = PjWj with (3j > 0, wj G G and uij C z 
for all j. As z is not C-minimal, this representation is not trivial. 

Note that in z = Y2, a i v ii z = T2Pj w j> an d in any intermediate representation z = YIPkUk, we have 
that the summands Vi, Wj, Uk have a strictly smaller norm on the first d components as z. In fact, the 
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same property holds for the sum of two summands that are iteratively needed for the rewriting steps. 
Thus, their corresponding orbits also have a strictly smaller norm than || orbs A (z)||- We conclude that 
at the time reps A (orbs A (z)) was chosen from G sym and then added to G sym , representatives for all 
orbits orbs A (ifj) were already in G sym . Thus, normalForm(reps A (orbs A (z)), G sym ) must have returned 
in contradiction to the assumption that reps A (orbs A (z)) was added to G sym . 

We conclude G\Q(A) = and hence G = £(A). □ 



5 Computational Experiments 

In this section we report on computational experience with a few examples that have symmetry. These 
problem deal all with 3- way tables, that is, with k\ x It2 X k% tables of unit cubes, each containing an 
integer number. For each problem, the lattice A is the set of lattice points in the kernel of the matrix 
that encodes the conditions that the sums along each one-dimensional row parallel to a coordinate 
axis are 0. Example 11.11 presented the case of 3 x 3 tables. 

The following running times demonstrate that, as expected, our symmetric algorithm heavily speeds 
up the computation. (Times are given in seconds on a Sun UltraSparc III+ with 1.05 GHz.) 



Problem 


\Sa\ 


Size of Graver basis 


|G sym | 


Algorithm 13. 11 


Algorithm 14.11 


3x3x3 


1,296 


795 


7 


2 


1 


3x3x4 


1,728 


19,722 


27 


1,176 


9 


3x3x5 


8,640 


263,610 


61 


560,517 


526 


3x4x4 


6,912 


4,617,444 


784 




260,590 
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